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ABSTRACT 

Galaxy-galaxy lensing is a powerful probe of the relation between galaxies and dark 
matter halos, but its theoretical interpretation requires a careful modeling of various 
contributions, such as the contribution from central and satellite galaxies. For this 
purpose, a phenomenological approach based on the halo model has been developed, 
allowing for fast exploration of the parameter space of models. In this paper, we inves- 
tigate the ability of the halo model to extract information from the g-g weak lensing 
signal by comparing it to high-resolution dissipationless simulations that resolve sub- 
halos. We find that the halo model reliably determines parameters such as the host 
halo mass of central galaxies, the fraction of galaxies that are satellites, and their 
radial distribution inside larger halos. If there is a significant scatter present in the 
central galaxy host halo mass distribution, then the mean and median mass of that 
distribution can differ significantly from one another, and the halo model mass de- 
termination lies between the two. This result suggests that when analyzing the data, 
galaxy subsamples with a narrow central galaxy halo mass distribution, such as those 
based on stellar mass, should be chosen for a simpler interpretation of the results. 

Key words: galaxies: halos - methods: analytical 



1 INTRODUCTION 



Understanding the connection between the spatial distribu- 
tion of galaxies and dark matter (DM) is one of the most 
important problems in modern cosmology. From the per- 
spective of fundamental physics, this connection needs to be 
understood if we wish to use galaxy clustering to determine 
dark matter correlations, a potentially powerful discrimina- 
tor between cosmological models. This connection can also 
be a powerful test of the tenets of the standard cosmolog- 
ical paradigm, such as the collisionless nature of cold dark 
matter. From the astrophysics perspective, this connection 
is an essential ingredient in the physics of galaxy formation. 
Current models based on pure N-body simulations, semi- 
analytic models or hydrodynamic simulations can accomo- 
date many observational aspects of this connection, but sev- 
eral ingredients of these models remain uncertain and need 
to be inserted ad-hoc, so it is not clear how much of the suc- 
cess is a result of the allowed freedom within these models. 
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It is important to look for new ways to test this connec- 
tion, to confirm and improve the existing models of struc- 
ture formation, and to distinguish astrophysical effects from 
cosmological ones. 

One of the probes of the galaxy-DM connection that re- 
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tortions of background galaxies around foreground galaxies, 
allowing direct measurement of the galaxy-dark matter cor- 
relation function around galaxies. The individual distortions 
are small (of order 0.1%), but by averaging over all fore- 
ground galaxies within a given subsample, we obtain high 
signal to noise in the shear as a function of angular separa- 
tion from the galaxy. If we know the galaxy redshifts, the 
shear signal can be related to the projected mass density as 
a function of proper distance from the galaxy. This allows 
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us to determine statistically the dark matter distribution 
around any given galaxy sample. 

In recent years, the progress on the observational side 
of g-g lensing has been remarkable. In the latest Sloan 
Digital Sky Surv ey (SDSS) analyses JSheldon et alJ 12004 
ISeliak et al.ll2005l) . 20-30 sigma detections of the signal as a 
function of physical separation have been obtained. Similarly 
high S/N detections have also been observe d as a function 
of an gular separation with other surveys (Hoekstra et al. 
12004) . This increased statistical power has been accompa- 
nied by a more careful investigation of systematic errors, 
such as calibration biases and intrinsic correlations, which 
for the SDSS are currently around 10% and the refore already 
dominate the error budget jHirata et al J2004 . As the qual- 
ity and quantity of the data and its analysis improves, the 
reliability of theoretical interpretation must be improved as 
well. The goal of the present paper is to compare the various 
theoretical analyses among themselves, and to discuss how 
they should be applied to the data. 

Theoretical analysis of g-g lensing falls into two cate- 
gories. The first ap proach is a direct comparison of simu- 
lations to the data iGuzik fc Seliakll200li lYang et aljEoo3 
I Weinberg et alJl2004 iTasitsiomi et aljl2004) . This approach 
is direct, but rather expensive, since the process of galaxy 
formation is not understood sufficiently well to result in 
unique predictions for a given cosmological model, nor is the 
cosmological model itself determined yet. Moreover, current 
simulations still suffer from a limited dynamical range, in 
the sense that they require a high mass and force resolu- 
tion to resolve individual galaxies and their associated dark 
matter halos, while at the same time they must also have 
sufficiently large volume to simulate a representative region 
of the universe. Several simulations of varying box size are 
needed to cover the whole observational range in luminosity 
and scale. As a result, many different simulations would be 
needed to explore the whole range of parameter space; at 
present, only a handful of simulations have been used for 
this application. 

The second approach is to u se the halo model (for a re- 
view, see ICoorav fc Shethll2002l) to describe the connection 
between galaxies and dark matter IIGuzik fc Se liak 2002, 
hereinafter GS02). This approach is more phenomenologi- 
cal, and for a given galaxy class leads to determination of 
quantities such as the virial mass distribution and the frac- 
tion of these galaxies that are satellites, which are useful 
quantities for constraining the galaxy formation models and 
cosmological models. While the information extracted from 
the data using the halo model may be all that is needed 
to quantify the galaxy-dark matter connection, halo models 
are reliable only to the extent that they are able to repro- 
duce the simulations, so they must be tested and calibrated 
on the simulations before applications to the real data are 
reliable. 

In this paper, we attempt to underst and how wel l the 
phenomenological halo model developed in ISeliakU2000D and 
studied in GS02 can reproduce the simulations. We are inter- 
ested in its ability to reproduce the g-g lensing signal, which 
in turn can be used to determine the halo mass probabil- 
ity distribution, the radial profile of satellite galaxies inside 
larger halos, and other quantities of interest. 



2 SIMULATIONS 

We use simulations described in ITasitsiomi et all ]2004h . 
where the reader can find further details. We assume the 
concordance flat ACDM model: fi m = 1 — JIa = 0.3, h = 0.7, 
where Q m and Q,a are the present-day matter and vacuum 
densities, and h is the dimensionless Hubble constant de- 
fined as Ho = WOh km s _1 Mpc -1 . The power spectrum 
normalization is ag = 0.9. The mode l is consistent with 
recent observational constraints (e.g.. Ispergel et alJ 1200 j : 
iTegmark. et alJl2004) . The effects of the power spectrum 
no rmalization, box size a nd cosmic variance were studied 
in ITasitsiomi et all i2004) using a range of simulation box 
sizes. Given current constraints, it is impossible to achieve 
the desired level of precision across the entire range of lu- 
minosities probed by the observations. Small box sizes can 
achieve sufficient mass resolution to resolve very small ha- 
los, but the sampling variance due to large scale fluctua- 
tions is large, which results in large errors on predictions 
at scales above a few hundred kpc. To reduce sampling 
errors, we wil l choose the largest bo x available from the 
simulations in ITasitsiomi et al] ll2004) . a 120 /i -1 Mpc box 
with 512 3 particles. The particle mass of this simulation is 
m p — 1.07 x 10 9 /i _1 Mq, so only halos with virial mass above 
a few times 1O 11 /i" 1 M are resolved. For this reason, we will 
restrict the analysis to galaxies brighter than M r = —19. 

The simulations were run usin g the Adaptive Re- 
finement Tree A^-body code (ART; IKravtsov et al] Il997l : 
iKravtsov fc Klvpinll999h . The ART code reaches high force 
resolution by refining all high-density regions with an au- 
tomated refinement algorithm. The criterion for refinement 
is the mass of particles per cell. The initial grid is 1024 3 
and the refinement criterion is level- and time-dependent. 
At the early stages of evolution (a < 0.65) the thresholds 
are set to 2, 3, and 4 particle masses for the zeroth, first, 
and second and higher levels, respectively. At low redshifts, 
a > 0.65, the thresholds for these refinement levels are set 
to 6, 5, and 5 particle masses. The lower thresholds at high 
redshifts are set to ensure that collapse of small-mass halos 
is followed with higher resolution. The maximum achieved 
level of refinement is L max = 8. As a function of redshift 
the maximum level of refinement is equal to L max = 6 for 
5 < z < 7, L = 7 for 1 < z < 5, L 

max 

> 8 for z < 1. The 
peak formal resolution is /i pea k < 1.8/i _1 kpc (physical). 

A variant of the Bound Den sity Maxima halo find- 
ing algorithm llKlypin et al.l Il999l is used to identify ha- 
los and the subhalos within them. The details of the algo- 
rithm and parameters used in the halo finder can be found in 
IKravtsov et al] i2004 . The main steps of the algorithm are 
the identification of local density peaks (potential halo cen- 
ters) and analysis of the density distribution and velocities of 
the surrounding particles to test whether a given peak corre- 
sponds to a gravitationally bound clump. More specifically, 
density, circular velocity, and velocity dispersion profiles are 
constructed around each potential halo center. The unbound 
particles ar e then remov e d iter atively using the procedure 
outlined in iKlypin et al] il999h . The final profiles are con- 
structed using only bound particles. We use these profiles 
to calculate properties of halos, such as the circular velocity 
profile V C i IC (r) = \J GM(< r)/r, and compute the maximum 
circular velocity V max . 

In this study, we distinguish between host halos with 
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centers that do not lie within any larger virialized system, 
and subhalos (or satellites) located within the virial radii of 
larger systems. We associate the former with central galax- 
ies and the latter with non-central galaxies or satellites. To 
classify the halos, we calculate the formal boundary of each 
object as the radius corresponding to an enclosed overden- 
sity of 180 with respect to the mean density of the universe. 
Note that we do not consider the center of a host halo to be 
a subhalo. Thus, host halos may or may not contain any sub- 
halos with circular velocity above the threshold of a given 
sample. The host centers, however, are included in cluster- 
ing statistics because we assume that each host harbors a 
central galaxy at its center. Therefore, the total sample of 
galactic halos contains central and satellite galaxies. The 
former have the positions and maximum circular velocities 
of their host halos, while the latter have the positions and 
maximum circular velocities of subhalos. In a cluster, for 
example, the brightest central galaxy that typically resides 
near the center would be associated with the cluster host 
halo in our terminology. All other galaxies within the virial 
radius of the cluster would be considered "satellites" asso- 
ciated with subhalos. 

Our galaxy sample is created by assigning realistic 
SDSS luminosities and colors to dark matter halos. To con- 
struct mock galaxy catalogs for comparison with observa- 
tions, one must define selection criteria for particular halo 
properties to mimic the selection function of the observa- 
tional sample as closely as possible. Halo mass is often used 
to define halo catalogs; e.g., a catalog can be constructed 
by selecting all halos in a given mass range. However, the 
mass and radius are poorly defined for the satellite halos 
due to tidal stri pping, which alters a halo's mass and phys- 
ical extent (see lKlypin et af1ll999l) . Therefore, we use the 
maximum circular velocity, V max , as a proxy for the halo 
mass. For isolated halos, V max and the halo's virial mass are 
directly related. For subhalos, V max will experience secula r 
decrease but at a relatively slow rate jKravtsov et al]|2004l) . 

To mimic the observational selection function, r-band 
luminosities are assigned to the halos as follows. We match 
the cumulative velocity function n(> Knax) of the halos to 
the SDSS observed r -band cumulative luminosity function 
jBlanton. et alJ [2003 ) . Note that n(> Vmax) includes both 
isolated host halos and subhalos. We use the r-band data 
since that band was used for SDSS spectra selection, has a 
more reliable luminosity function measurement observation- 
ally, and is the focus of most SDSS analyses. We first obtain 
the average Vmax — M r relation by matching n(> Vmax) to 
n (< Mr). This was the same method used to assign luminosi- 
ties to subhalos in lKravtsov et alJ 11200411 . in which galaxy 
clustering properties were reproduced remarkably well. The 
mean redshift of the lens galaxies is 0.1, so we use halo cat- 
alogs at this redshift. 

One may also introduce scatter in the relation between 
Knax and M r , which here we have chosen as a Gaussian in 
M r at fixed Vmax, as described in detail in iTasitsiomi et all 
(2004). The amplitude and dependence of the scatter on 
galaxy luminosity or halo mass are very uncertain, so here 
we will simply explore this example to investigate its conse- 
quences. The value of scatter is meant to be realistic for 
the current data, which are a mixture of early and late 
type galaxies and for which intrinsic reddening corrections 
have not been applied. The Gaussian scatter is introduced 



in a fashion that keeps the luminosity function constant, 
which results in somewhat lower scatter for higher luminos- 
ity galaxies and vice versa. Below we present extensive com- 
parisons of the halo model and the actual simulation results 
with and without scatter. 



3 OVERVIEW OF THE HALO MODEL 
3.1 Dark matter halos 

In current cosmological models, structure grows hierarchi- 
cally from small, initially Gaussian fluctuations. Once the 
fluctuations go nonlinear, they collapse into virialized halos. 
The spatial density of halos as a function of their mass M is 
specified by the halo mass function dn/dM, which in general 
is a function of redshift z. It can be written as 



^-dM = &f(u)du, 
dM M 



(1) 



where p is the mean matter density of the universe. We in- 
troduced the function f{y), which can be expressed in units 
in which it has a theoretically universal form independent 
of the power spectrum or redshift if written as a function of 
peak height 



(2) 



<r(M) 



Here S c = 1.686 is the linear overdensity at which a spherical 
perturbation collapses at redshift z, and a(M) is the rms 
fluctuation in spheres that contain on average mass M at 
an initial time, extrapolated using linear theory to z. 

The first analytic model for th e mass function was pro- 
posed bv lPress fc Schechterl 1)1974?) . While it correctly pre- 
dicts the abundance of massive halos, it overpredicts the 
abundance of halos around and below the nonlinear mass 
scale M n i(z) (defined below). An imp roved version has been 
proposed bv lSheth fc Tormenl <1999» . 



vf{v) = A(l + i/"*) 



P\ '1/2 



V/2 



(3) 



where v' = a v, with a = 0.707 and p = 0.3 as the best 
fit values (the IPress fc Schechterl lll974ll mass function cor- 
responds to a = 1, p = 0). Furt her modificat i ons to this 
expression have been suggested in lYahaei et alJ (12004ft . but 
the effects are very small and we will ignore them here. The 
constant A is determined by mass conservation, requiring 
that the integral over the mass function times the mass gives 
the mean density or, equivalently, J f(y)dv = 1. It has been 
shown that the form in equation |3] can be derived analyti- 
cally within the fra mework of the ellipsoidal collapse model 
JSheth et alJl200ll) . For a given peak height, we can also 
compute the bias 



b(y) 



1 



2p 



8 C 5 C (1 + v' p ) ' 



(4) 



where for the purpose of computing the bias we use a = 0.73 
and p = 0.15 rather than the values from the mass function, 
in ord er to best 
(2004). As shown in! 



ler than the values from t he mass function, 
match th e results from ISelj^k fc Warrenl 
in F ig. 2 o f lWarren. et alJ ll2005h . the origi- 
nal lSheth fc Tormenl i 19991) mass function and this modified 
one give nearly identical results in the most relevant range of 
halo masses, lO^-lO 14 Ii~ 1 Mq. In support of this assertion, 
we note that when we compare the predicted halo model 
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lensing signal computed with the two mass functions, the 
difference is less than 1 per cent. 

The halo mass is defined in terms of the linking length 
parameter of the friends-of-friend s (FoF) algorithm , whic h 
is 0.2 for the simulations used in ISheth fc Tormenl Jl99St) . 
This definition roughly corresponds to spheric al overdensity 
halos of 180 time s the background density djenkins et alJ 
l200ll:IWhitell200lD . For the range of masses of interest here 
and their corresponding halo concentrations, it is about 30- 
35% larger than the mass defined as the mass within the 
radius where the density is 200 times the critical density. 
Since we use fl m = 0.3 when computing the virial masses, 
they are defined as the mass within the radius in which the 
mean density is 54 times the critical density. We will also 
define the concentration parameter relative to this radius. 



3.2 Halo-galaxy connection and galaxy-galaxy 
lensing 

In the current paradigm of structure formation, all galaxies 
form inside dark matter halos. While this is generally ac- 
cepted, we also know that the relation between the two is 
not one-to-one, and galaxies of the same luminosity can be 
found in halos of different masses. For example, a typical 
galaxy like the Milky Way may be found at the center of a 
low mass halo with a typical size 200 kpc, it may be part of a 
small group with typical size 500 kpc or it may be a satellite 
in a cluster with a typical size of 1-2 Mpc. We wish therefore 
to determine the probability AP for a galaxy in this sample 
to be in a halo of mass M ± AM/2. To describe this prob- 
ability, we will use the conditional halo mass probability 
distribution dP/dM = p{M;L) (GS02). With g-g lensing, 
one can in principle determine the full halo mass function, 
since small halos contribute only at small scales, while large 
halos such as clusters also generate signal at larger scales. 
The fact that g-g lensing measures the signal over a wide 
range of scales facilitates the determination of the full halo 
mass function for a given subsample. In practice, given the 
measurement errors, a model must be adopted to extract 
this information from the data. 

Galaxy-galaxy lensing measures the tangential shear 
distortions in the shapes of background galaxies induced by 
the mass distribution around foreground galaxies. Because 
the shear distortions 74 are very small, in our case 10 -3 , 
while the typical galaxy shape noise is 0.3, we must average 
over many foreground-background pairs to extract the sig- 
nal. The result is a measurement of the shear-galaxy cross- 
correlation as a function of relative foreground-background 
separation on the sky. We will assume that the redshift of 
the foreground galaxy is known, so one can express the rel- 
ative separation in terms of transverse comoving scale R. 
If, in addition, the redshift distribution of the background 
galaxies, or their actual redshifts, are known, then one can 
relate the shear distortion -y t to AE(i?) = E(< R) - E(i?), 
where E(J?) is the surface mass density at the transverse 
separation R and E(< R) its mean within R, via 



7* 



E C rit 



Here 



Ecrit — 



rs 



47rG (1 + z L )r L r L s ' 



(5) 



(6) 



where tl and rs are the comoving distances to the lens 
and source, respectively, tls is the comoving distance be- 
tween the two and zl is the redshift of the lens. If only 
the probability distribution for source redshifts is known, 
then this expression needs to be integrated over. In princi- 
ple, the relation between comoving distance and measured 
redshift depends on cosmology, but since we are dealing with 
low redshift objects, varying cosmology within the allowed 
range makes little difference. We will assume a cosmology 
with £l m — 0.3 and J1a = 0.7, and we work in comoving 
coordinates throughout the paper. 

We will first overview the halo model formalism of 
GS02, beginning the discussion with a simplified description. 
Let us assume that a given halo of mass M produces an av- 
erage lensing profile AE(J?, M). This can be obtained via 
line of sight integration over the dark matter profile, which 
in th is paper is modeled as a NFW profile jNavarro et alJ 
1996) 



p(r) = 



(r/r s )(l + r/r s y 



(7) 



This model assumes that the profile shape is universal in 
units of the scale radius r s , while its characteristic density 
p s or concentration Cd m = r v /r s may depend on the halo 
mass, which he re will be modeled as feullock et al.l l200ll : 
lEke et al.ll200l 



f-'dn 



10 



A I 



M nl (z) 



(8) 



M n i(z), which depends on the cosmology, is defined such 
that the rms linear density fluctuation, extrapolated to red- 
shift z, within a sphere containing mass M n i is equal to 5 C . 
Since most of the signal is at R > 50 — 100 /i~ 1 kpc. bary - 
onic effects can be neglected (see, e.g.. lGnedin et ai"ll2004l) . 
dark matter profiles are well determined from simulations, 
and concentration or the choice of the halo profile does not 
play a major role. The average g-g lensing signal for a galaxy 
with luminosity L is 



(AE)(R;L) = p(M;L)AT,(R,M)dM. 



(9) 



Because of noise in the g-g lensing signal, we cannot invert 
the relation to obtain the conditional mass probability dis- 
tribution with arbitrary precision, but instead assume some 
functional form for it and fit for its parameters. 

Our simplified description so far ignores the fact that 
there are two distinct galaxy types that need to be modeled 
separately. The first type are the galaxies that formed at 
the centers of dark matter halos, such as the so-called field 
galaxies or cD galaxies in the cluster centers. The second 
type are the non-central galaxies, or satellites. We know that 
a galaxy of a given luminosity can be of either type, so fol- 
lowing GS02, we split p(M; L) into two parts, p c and p NG , 
representing respectively central and non-central galaxies, 
with the fraction of non-central galaxies in each luminosity 
bin Li given by a free parameter on: 



V{M-U) 



1 



a l )p c {M;L,)+a l p NC {M;L l ) . (10) 



on is, by definition, the satellite fraction determined for the 
given lens sample as a whole. 

For the central galaxy population we assume that the 
relation between the halo mass and galaxy luminosity is 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy-galaxy lensing 5 



tight and we model this component with a delta-function, 

p c (M;Li)dM = S D (M - M ,i)dM , (11) 

where Mo t i is the typical halo mass of the ith luminosity bin. 
In reality, this component should have some width both be- 
cause of intrinsic scatter in the M — L relation, scatter due 
to inclination and dust, and because we work with luminos- 
ity bins of finite width. Explicit tests have shown that the 
results are only weakly affected even if the scatter is a fac- 
tor of 2-3 in mass (GS02). If the scatter is larger, then the 
effects can be larger; we will use simulations to investigate 
this effect. 

The non-central galaxies are different in that they have 
presumably formed in smaller halos which then merged into 
larger ones. It is thus reasonable to assume that their lumi- 
nosity is not related to the final (host) halo mass. Instead we 
assume a relation between the number of these non-central 
galaxies and the halo mass: the larger the halo, the more 
satellites of a given luminosity one expects to find in it. We 
assume this relation is a power law, {N)(M; L) oc M E , above 
some minimal halo mass Af m i n , which should be larger than 
the halo mass of the central galaxy component above, since 
we are assuming that there is already another galaxy at the 
halo center. Below this cutoff, the number of galaxies quickly 
goes to zero. These assumptions imply 

P NC (M; LiidM oc F(M)M H —dM . (12) 

dM 

GS02 used F(M) = e H (Af-M minii ) where M minji = 3M 0yi ; 
here, we will use a slightly more realistic functional form 
that is smoother than the step function and matches the 
simulations better. This is described in more detail in the 
next section. We have verified that the two expressions do 
not yield significantl y different results. Semi- analytic models 
of galaxy formation (iKauffmann et alJll99Sl GS 02), subha- 
los in N-body simulations IIKravtsov et alJl2004l). and direct 
obser vational measurements jLrnetall]20f5[Zehavi. et alJ 
l2004f) agree with this model, and predict that for most galax- 
ies, e ~ 1 and a ~ 0.2. In fact, when attempting to fit for 
both e and a, we find that these two parameters are almost 
completely degenerate. We therefore assume fixed e = 1, and 
do not fit for it throughout this study. 

For the non-central component, the weak lensing profile 
AE(_R, M) is a convolution of the halo profile with the radial 
distribution of the galaxies. Since we are explicitly excluding 
the central galaxies, the non-central galaxy component of 
the g-g lensing signal does not peak at the center. Instead, 
for a given halo mass, it is small at small radii, peaks at a 
fraction of the virial radius, and drops off at large radii. 

We assume that the radial distribution of galaxies is 
proportional to the dark matter profile, c g = acd m - Obser- 
vationally, there is not much evidence for any departures 
from a — 1 , but the constraints are weak and only exist 
for clusters llCarlberg et al]ll997t) . As a result, a could dif- 
fer sig nificantl y from un ity, as discussed in more detail in 
iNaeai fc Kravtsovl i2005l) . We will show that g-g lensing may 
be a useful way to measure it observationally for groups as 
well. 

On scales below 2/i _1 Mpc, the one-halo or Poisson term 
dominates the g-g lensing signal. As shown in GS02, the lens- 
ing signal from neighboring halos (called the "halo-halo" or 
hh term) can be neglected for R < 2/i _1 Mpc, while on very 



large scales this term dominates and follows linear theory. 
For the sake of consistency between results on large and 
small scales, we have included it here. Since the NFW pro- 
file is not abruptly truncated at the virial radius in simu- 
lations, we have investigated the dependence of the results 
on the maximum radius to which the halo profile is inte- 
grated in, e.g., equation 1141 We define this radius in units of 
the virial radius, so x\ = 1 means that no mass beyond the 
virial radius is taken into account. (This is only true for the 
non-central term; when computing the Poisson term due to 
the central halo, the profile is not truncated at xir V i T when 
projecting to get E(ii) from p(r).) Since we do not assume 
that the halo mass function accounts for all the mass in the 
universe, there is no concern of double-counting mass when 
using x\ > 1. The halo-halo correlations also depend on the 
large scale bias, which for a given halo mass distribution is 
given by equation^] Unfortunately, the small box size of the 
simulations implies that the sampling variance is important 
on large scales, so the errors on large scales may be underes- 
timated. They are also highly correlated, and our fits do not 
take the correlation into account. Nonetheless, we perform 
fits out to 13 h~ Mpc, keeping in mind that the errorbars 
are actually larger due to the correlations, and that the bias 
is somewhat degenerate with the choice of xi . Based on our 
investigations, we see no reason to include mass beyond the 
virial radius, and so use xi = 1 throughout this work, but 
larger simulations will be necessary to investigate the be- 
havior on large scales in more depth. 

In addition to correlations of satellite galaxies with the 
host halo dark matter, there are correlations between a satel- 
lite galaxy and dark matter particles of its own (sub)halo. 
Thus, the remaining uncertainty in modeling satellites is 
how much of the dark matter around satellite galaxies re- 
mains attached to them. Analyses of cosmological simula- 
tions show that the typical fraction of mass bound to subha- 
los is « 0.1 — 0.2, roughly independent of the host halo mass 
but stron gly dependent on the distance from the host cen- 
ter ( e.g., iGhigna et ai]ll99St iGao et aLlliool IWeller et all 
2004). Since the fraction of subhalos, a, is typically low 
(a ~ 0.2), the correction due to the subhalo bound mass 
is small. The fraction of dark matter remaining is likely to 
depend on the details of satellite history, such as the point 
of closest approach to the halo center, where the tidal strip- 
ping is strongest. It is also likely to be somewhat correlated 
with the instantaneous position of the satellite, but this cor- 
relation is not perfect, given the highly elliptical orb its of 
satellites observed in simulations dGhigna et al.ll2000t) . De- 
termining this position in all but the most massive clusters 
(where the center is often traced by either X-ray emission or 
a cD galaxy) is difficult. For this reason, we will work with 
an average quantity rather than some function of radial dis- 
tance from the halo center. We assume the dark matter was 
tidally stripped in the outer parts of the halo, but remains 
unmodified in the inner parts of the satellite halo. Effec- 
tively this means that each non-central galaxy also has a 
central contribution, which we model in the same way as 
for the central galaxies (i.e., as a halo with mass Mo,% be- 
fore stripping) out to a fraction of virial radius, but totally 
stripped beyond that, yielding g-g lensing signal AE oc R~ 2 . 
Based on fits (figure Ujjl we will choose the truncation radius 
at 0.4r v i r , which is equivalent to having 50% of the mass 
stripped, consistent with the average mass loss for subha- 
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los observed in cosmolo gical simulations llGao et al]l2004t 
iNaeai fc Kravtsovll2005h . We note, however, that given ex- 
isting measurement errors, the g-g lensing signal is only 
weakly sensitive to the exact value. Note that this trun- 
cation is responsible for the kink in the stripped satellite 
profile in the noncentral signal in figure Q a more realistic 
approach might involve smoothing out this feature. 



3.3 Signal computation 

Using all of the above ingredients, we compute the predicted 
signal AE(_R, M) for a given cosmological model and central 
halo mass Mo, central galaxy bias b(Mo), and c g . Here we 
briefly describe the signal computation for the various con- 
tributions. Note that when we compute signal in comoving 
coordinates, the characteristic redshift zo is only used when 
computing the growth factor b(zo) for use in normalizing the 
linear power spectrum, and when determining M n \(z). 

First, there is the 1-halo (Poisson) term for the central 
galaxies, AEf . This term can simply be computed by finding 
E(i?) and then AE(ii) for a NFW profile given M and 
Cdm which is determined from equation [H] When combining 
the signal from central and satellite galaxies, this term is 
multiplied by 1 — a. There is also a similar term for the 
satellites, AE^ t , which assumes a stripped density profile. 
This term is the same as that for the central galaxies out to 
0Ar v i r (our chosen truncation radius), then is proportional 
to R~ 2 beyond that scale. When combining signal for central 
and satellite galaxies, this term is multiplied by a. 

For the remaining terms in the signal, AE was com- 
puted by finding the galaxy-dark mat ter cros s powe r spec- 
trum (this approach is described in ISeliakl (2000)). The 
power spectrum can then be Fourier transformed to ob- 
tain the galaxy-dark matter cross-correlation function £g,dm, 
which yields E(J?) via an integration over the comoving sep- 
aration 



x 



2x1/2 



]dx 



(13) 



(where we have dropped an unobservable constant term) and 
AE(ii) after another integration. The power spectrum can 
be related to the Fourier transform of the halo profile via 
y(k, M), where 



y(k,M) = 



pjk, M) 
M 



1 

M 



2 , , . sin (kr) 

4vrr ■ drptr, M) ± — '- 

kr 



(14) 

We can define j/ g , the Fourier transform of the radial profile 
of galaxies using a NFW profile with concentration c 9 , and 
ydm, the Fourier transform of the halo dark matter profile 
with concentration Cd m - 

This procedure was used to get the noncentral Poisson 



Pnc(k) = 



f{u)du(N{M))y dm {k, M)y s (k, M) 



(15) 

and describes the correlations between satellite galaxies and 
their host halo dark matter. Because this is a noncentral 
term, it must be multiplied by a when computing its con- 
tribution to AE. 

This procedure was also used to compute the central 
and noncentral halo-halo terms, AEj^t and AE„c , which 



dominate on large scales. The power spectrum for the central 
halo-halo term (describing correlations between the central 
galaxy of one host halo and the dark matter of another host 
halo), is 

PcLik) = b(M )P Un (k) J f{v)dub(u)y dm (k,M) (16) 

where the integration over v includes all values. Note that 
dark matter must on average not be biased, and therefore 
the last term for a general b(v) must reduce to 1 on large 
scales where jfc = 1. For the bias expression used in this 
paper, equation^] the average bias is 1.06 rather than 1, but 
in order to reproduce the relation -P C cnt(&) = b(Mo)Pu n (k) 
on large scales, we have artificially normalized the integral 
to 1 on large scales. This central halo-halo power spectrum 
must be multiplied by (1 — a) when combining the signal for 
a mixture of central and noncentral galaxies. 

For the remaining noncentral halo-halo term, we have 

P^(k) = P lm (k) J f(v)dvb(v)y Am (k,M)x 



£ I f{»)dub{v) 



(N(M)) 
M 



y g (k,M), (17) 



where the first integral, as above, must be 1 on large scales, 
and the second term can be normalized via 



f(v)du 



(N(M)) 
M 



and (once normalized) equals the average bias (b) for satel- 
lite galaxies of this mass. Consequently, Pnc(k) reduces to 
(b)Piin (k) on large scales, and must be multiplied by a when 
combining signal for central and satellite galaxies. When the 
central and noncentral signal are combined, we then have on 
large scales an effective bias 6 e ff = (1 — a)b(Mo) + a{b). 

To illustrate the effects of the various contributions to 
the signal discussed in this section, figure Q shows the pre- 
dicted central (top), noncentral (middle) and total (bot- 
tom) signal for Q, m = 0.3, erg = 0.9, b(M ) = 0.8, and 
Mo = 5 x 10 12 h" 1 Mq. The central signal is the sum of 
two contributions: 



AE ccn t — AE ccnt + AE ccnt 



(18) 



and the satellite signal is the sum of three contributions: 



AE sat = AEf at + AE„ C + AEnc- 



(19) 



For the combined signal we used satellite fraction a = 0.2. 



4 HALO MODEL VERSUS SIMULATIONS 

As described in the introduction, the main motivation for 
the comparison between the halo model and simulations is 
that the halo model makes several approximations, so it 
must be verified and calibrated on realistic simulations. We 
may, of course, find that the simulations reproduce the data 
perfectly, in which case one could argue that the halo model 
itself is unnecessary. However, there may be different cos- 
mological models or different treatments of how to populate 
halos and subhalos with galaxies that all agree equally well 
with the data. Rather than identify them with many expen- 
sive numerical simulations, one can use the halo model to 
do the same. 
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Figure 1. Halo model predictions for AE(r) for a central halo 
mass of Mo = 5 X 10 12 fe _1 Mq. The top panel shows the cen- 
tral signal with the 1-halo and halo-halo terms shown separately, 
assuming b(Mo) = 0.8. The middle panel shows the noncentral 
signal, with the stripped satellite central profile, the 1-halo non- 
central term, and the halo-halo noncentral term shown separately 
as labeled. The bottom panel shows the combined signal, with 
satellite fraction a = 0.2. 



As discussed in §2, the simulations used here cannot 
reliably resolve low mass halos, so we restrict our analysis 
to 3 luminosity bins, each 1 magnitude wide in r-band, L3 
[-19, -201 L4 [-20,-2 1 1 and L5 [-21, -22], using the no- 
tation from lSeliak et all (2005). We do not use the brightest 
bin (L6 [—22,-23]), since most of these galaxies reside in 
clusters, and our simulation box is too small to properly 
sample these. 

To test how well g-g lensing analysis can extract various 
halo model parameters, we take the full g-g lensing signal 
from simulations and fit for 4 parameters: Mo, b(Mo), satel- 
lite fraction a, and radial distribution of satellites c 9 . When 
performing the fits to the central galaxy signal alone, we 
also explored the sensitivity to the dark matter halo con- 
centration parameter Cd m . We assign errors to the radial 
poi nts using the actu al errors that we find in SDSS analy- 
sis (ISeliak et all2005l) . As the survey continues to take data, 
the error amplitude will decrease, but the relative errors will 
remain the same, so one can simply rescale the errors by the 
appropriate a mount. Note, howeve r, that a deeper survey, 
such as RCS IIHoekstra et al.ll2004F) . can provide better in- 
formation on g-g lensing at small radii, which would be par- 
ticularly useful for determination of the halo concentration 
parameter Cd m . 

Figure [5] shows the lensing signal AE(i?) as a function 
of transverse separation R from the simulations, together 
with the halo model signal that fits the simulations best. 
We show the results separately for all luminosity bins with 
and without scatter. The fits are done using SDSS error bars 



JSeliak et alJl2005l) reduced by a factor of 10, with the error 
from simulations added in quadrature. 

To obtain an estimate of the minimum reliable scale in 
simulations, we performed a convergence test with respect 
to AE, using the low and high resol ution (8 times more par - 
ticles) simulations of cluster CL2 of lTasitsiomi et all ll2004i) . 
We found that a reliable minimum scale is ~50 /i _1 kpc and 
hence for the fits we only use simulation data on scales larger 
than that. In fact, resolution can affect both the dark matter 
density profiles, making them somewhat shallower than they 
actually are, and the number of subhalos in the innermost 
host halo regions. 

The most striking feature of the figure is how well the 
halo model reproduces the simulations with 3 free parame- 
ters (central halo mass, satellite fraction and satellite radial 
distribution, assuming that bias is determined by the halo 
mass, which is confirmed when we compare the fitted values 
to those predicted). As discussed, the small disagreement on 
small scales is actually caused by finite resolution in simu- 
lations, which leads to an enhancement of AE. 

Figure shows the signal from simulations for a com- 
bined sample of central and satellite galaxies, such as would 
be observed in a typical luminosity-selected sample with no 
additional selections applied. We have also done the fits sep- 
arately for central and noncentral components. If one takes 
galaxies in dense regions, then the fraction of satellites is 
likely to be increased, while if one takes the sample of galax- 
ies in low density regions, the satellite fraction is likely to 
be decreased. Similar effect can be achieved by color selec- 
tion; for example, when taking blue galaxies, the satellite 
fraction is likely to be decreased. These two sub-samples 
therefore provide the extreme cases of such selection, and 
are useful as an indication of what additional information 
one may be able to extract from the g-g lensing data given 
such selections. 

4.1 Central galaxies 

There are several aspects of the halo model that we wish to 
test. We begin with the halo mass distribution for central 
galaxies. For each of the magnitude bins we have a distribu- 
tion of central halo masses from the simulations, as shown 
in figure 13 As is clear from this figure, the halo mass dis- 
tributions are reasonably narrow in the simulations without 
scatter, while in the case with scatter they are significantly 
broader, with the width and assymetry of the distribution 
increasing for higher luminosity bins. This effect is also evi- 
dent from the difference between the mean and the median 
of the mass distributions given in table which reaches a 
factor of 5 in the brightest bin with scatter. The width is 
largest for the brightest bin because the galaxies are on the 
exponential slope of the luminosity function, while many of 
the associated halos are on the exponential tail of the halo 
mass function. In light of the large width, it is difficult to 
speak of a single mass associated with a given luminosity 
bin, yet we attempt to determine a single mass using the 
fitting rather than the full shapes of the curves in Fig. |H] so 
we must investigate further to understand the consequence 
of this choice. The relevant questions associated with the 
halo model fitting are: What is the meaning of the best fit 
halo mass from the halo model fits? How much of an error is 
one making by using the best-fit masses if one is interested 
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r [Mph/h ] 

Figure 2. AE(r) for L3, L4, and L5 luminosity bins, bottom to 
top. The top panel is for simulations without scatter between 
Knax and L r , while the bottom panel is with scatter. Points 
are the simulation data with errors taken from observations and 
scaled down by a factor of 10, with simulation errors added in 
quadrature. The solid line is the halo model fit. Resolution effects 
become important below r~50 h~ J kpc. Thus, for the fits we used 
simulation data above this scale only; the disagreement for the in- 
ner few points is caused by insufficient resolution in simulations. 



in the mean (median) halo mass for a given sample? Can 
one identify from the lensing signal itself the width of the 
central halo mass distribution? 

To address the first two questions, we can compare the 
mean and median of the distribution to the best-fit masses 
using a NFW profile. We perform NFW fits in two sepa- 
rate ways. First, we fit the central galaxy component signal 
only to a NFW profile, using all of the information out to 
260ft -1 kpc. We chose this scale motivated by observational 
studi es where this scale i s used to define the aperture mass 
fe.g.. lMcKav. et all200lh . We fit simultaneously for concen- 
tration and virial mass. In the second case, we use the full 
combined signal, effectively cutting out the information from 
the central halos at large transverse separations, since the 
signal there is dominated by satellites sitting in more mas- 
sive halos, whose signal dominates at larger radii. In this 
case we assume concentration Cd m based on equation El and 
also fit simultaneously for b(Mo), c g and a. Both fits give 
similar results for virial masses (table 0. 

If the central galaxy halo mass distribution is narrow, 
the halo model is able to determine it quite accurately, as 
shown in tabled If the halo mass distribution is broad, then 
there is no typical mass, and in general the g-g lensing mass 
determination from NFW fits falls between the mean and 
the median mass of the halo population. If one is interested 
in the mean halo mass, then the NFW fits underestimate 
the mass by 10-30% at L3 and 10-55% at L4 and L5, so one 
must increase the best-fit NFW masses by up to a factor of 
2 or even more at the bright end when scatter is significant. 



10" 10 1Z 10 13 10" 




Figure 3. Halo mass distribution for central galaxies for L3 
(solid), L4 (dotted), and L5 (dashed) luminosity bins. The top 
panel is for simulations without scatter, the bottom is with scat- 
ter. Median and mean masses are given in table 1. 

If the median masses are of more interest, then one must de- 
crease NFW masses by similar amounts. At low luminosities, 
these corrections are small, but become increasingly impor- 
tant for brighter galaxies because the halo mass distribution 
is broader. Even more problematic is the fact that it is dif- 
ficult to assess the corrections from the fits itself, since the 
fits are good in both cases, but the corrections vary from 
essentially none in the case without scatter to over a factor 
of two in the case of scatter. To improve the reliability of 
the mass determination, one must therefore use samples of 
galaxies with as little scatter as possible. Below we discuss 
this point further. 

Despite the wide distribution of halo masses, the NFW 
fitting is able to extract the mean concentration parameters 
of the dark matter halos correctly. This is particularly true 
for the lower luminosity cases L3 and L4, where equation 
with M n i = 9 x 10 12 /i _1 Mq {Q, m = 0.3, a 8 = 0.9) gives 
Cdm = 12 — 14 for the mean virial mass in these bins, in good 
agreement with the fitted value Cdm ~ 13. In L5 without 
scatter, the fitted value, Cd m ~ 12, is somewhat higher than 
the predicted value Cdm ~ 10. It is possible that for this bin, 
the difference in best-fit masses between the two kinds of fits 
stems from the assumption of a lower concentration in the 
fit to the full signal, which would encourage a higher best-fit 
mass, as is observed. We note here, however, that the NFW 
dark matter profile is not a perfect fit to the simulations, so 
there are small differences in the NFW fits to 3-d data versus 
2-d AE data. Hence, these differences may be a reflection 
of the uncertainties in NFW parameterization rather than a 
failure of the halo model. 

Adding scatter leads to a modest decrease in the fit- 
ted concentration parameters. This effect is relatively small 
in the faintest bin and increases towards the brighter bins, 
where scatter has more effect. Thus, averaging over 1 mag- 
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Table 1. Mean and median halo masses for central galaxies in simulations with and without scatter for 3 luminosity bins. Also given are 
the mass and concentration from NFW profile fits to the g-g lensing signal for the central component only, and using the full (central and 
satellite) g-g lensing signal. For concentration C(j m , we only show the values for the fits to the central component; we assume equation 181 
for the fit to the full signal. 



Luminosity (M) centra j 


-^median, central 


A^NFWfit, central c NPWfit .central 


A%FWat,full Scatter in V m ax — L 


lO 12 /!" 1 M 


lO 12 ^- 1 M 


lO 12 /!" 1 Mq 


lO 12 ^ 1 Mq 



L3 [-19,-20] 0.506 0.460 0.476 12.8 0.42 ± 0.03 no 

L4 [-20,-21] 1.79 1.47 1.61 12.7 1.52 ±0.07 no 

L5 [-21,-22] 13.2 8.36 9.71 12.2 12.6 ± 0.3 no 

L3 [-19,-20] 0.745 0.390 0.579 11.6 0.47 ± 0.04 yes 

L4 [-20,-21] 2.91 1.10 1.86 10.9 1.42 ±0.09 yes 

L5 [-21,-22] 11.7 2.34 6.21 10.3 4.8 ± 0.2 yes 



Table 2. Fit results from a fit for central halo mass and bias b(Afo) to the signal for central galaxies, using information out to 13h 1 Mpc. 
Because neighbor ing bi ns are correlated on large scales, the bias errors are likely to be underestimated. The true bias values are from 
ISeliak fc WarrerJ <2004l) . 



Luminosity A^NFWfit, central ^(^ccntral) ^true Scatter in V m 



lO 12 ^- 1 Mr. 



L3 


[-19,-20] 


0.45 ± 0.02 


0.65 ±0.04 


0.68 


no 


L4 


[-20,-21] 


1.62 ±0.03 


0.68 ±0.04 


0.72 


no 


L5 


[-21,-22] 


11.8 ±0.2 


0.90 ±0.05 


1.00 


no 


L3 


[-19,-20] 


0.54 ±0.02 


0.64 ±0.04 


0.68 


yes 


L4 


[-20,-21] 


1.82 ±0.05 


0.73 ±0.04 


0.72 


yes 


L5 


[-21,-22] 


6.8 ±0.2 


1.00 ± 0.06 


0.93 


yes 



nitude bins in luminosity should allow one to extract the 
concentration parameter from g-g lensing signal with little 
or no bias for galaxies below L*. This is particularly the 
case if one can identify a galaxy population with a narrow 
distribution of halo masses. An example of this would be 
selecting elliptical galaxies that lie in a narrow strip on the 
fundamental plane which co rresponds to the dynamical mass 
iPadmanabhan et al.ll20(5l. Another example ar e galaxies 
selected by stellar mass iKauffmann. et al.1 12003). which is 
likely to be a more faithful tracer of dynamical mass than 
the luminosity itself. Careful selection of the lens population 
in one of these ways, via a reliable tracer of the dynamical 
mass, should help ensure that the best-fit central halo mass 
has a simple interpretation. 

If one assumes a given cosmological model, then the 
concentration parameter is fixed, and the difference between 
the predicted and observed concentration can be used to 
study the width of the halo mass distribution. We see that 
even in the cases of a very broad distribution, as in L5, 
the differences in the average concentration parameter are 
modest and it will be difficult to observe them directly in 
observations. Moreover, one can change concentration for a 
given halo mass by changing the cosmological parameters 
(primarily nonlinear mass), so the two effects are, to some 
extent, degenerate. 

Table [5] shows results for fits to the central signal us- 
ing all scales (out to 13h~ Mpc) to determine Mo and bias 
b(Mo). This determination of b is particularly useful since on 
large scales the central signal is completely negligible com- 
pared to the halo-halo term, so the bias is easily extracted 



from the signal amplitude. However, the errorbars are likely 
underestimated due to the correlations in the correlation 
function bins. As shown in table the bias values are rea- 
sonable for central h alos with those masses w hen compared 
against the results in ISeliak fc WarrerJ faOO^i . 

4.2 Satellites 

We turn next to the satellite signal. We begin with the dis- 
tribution of host halo masses for satellites in a given narrow 
luminosity range. These are very broad, as shown in figure 
f4] This is expected, since a satellite of a given luminosity can 
be part of a small group or a large cluster and there is little 
relation between the halo mass and satellite luminosity. The 
latter is particularly true for the case with scatter, where 
the differences between the distributions are very small and 
only appear at the low mass end of the host halo masses. 

Instead of discussing the actual mass distribution we 
can phrase it in terms of the average halo occupation model, 
(N(M)). To obtain the actual halo mass distribution of 
satellites, one has to multiply (N(M)) with the halo mass 
function (equation [3J|. The latter is well determined from 
simulations for masses in the range (10 12 — 10 14 )fc~ 1 M l ^, the 
range of interest here dSheth fc Tormenlll999t Ijenkins et alJ 
2001). The halo occupation (N(M)) as a function of halo 
mass M is shown in figure [S] for the 3 magnitude bins. We 
show simulation results without scatter and with scatter. 
We see that in all cases (N(M)) oc M (i.e., e — 1) is a good 
fit to the simulations above 3Mo, while below this mass the 
number of satellites declines more rapidly. The distributions 
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Figure 4. Halo mass distribution for hosts of satellite galaxies 
for L3 (solid), L4 (dotted), L5 (dashed) luminosity bins. The top 
panel shows the case where there is no scatter between V ma x and 
L r , while the bottom panel assumes scatter, as described in § 121 
Note that scatter makes the distributions nearly independent of 
luminosity. 
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Figure 5. Mean halo occupation number as a function of halo 
mass for satellite galaxies for the three luminosity bins, as la- 
belled. Solid triangles are without scatter, open squares are with 
scatter. The solid line is the analytic approximation used here for 
the halo model, the vertical line delineates the break at 3Mo as 
used in GS02. 



with and without scatter are very similar for L3 and L4, 
while for L5, the scatter increases the abundance of galaxies 
at the low mass end. We also show the simple step function 
model of GS02 and an improved model where the decline 
below 3Mq is more gradual, modeled here as (N(M)) oc M 2 
below 3Mo (the improved model was used for all fits). Note 
that the overall amplitude is not relevant in these relations, 
since we normalize it to the overall fraction of the galaxies 
that are satellites in any given sample. 

Examination of figure [5] on ~ 1 Mpc scale shows that 
scatter has little effect on the noncentral signal in L3 and 
L4, while for L5 the effect is much more significant: scatter 
increases the abundance of low mass halos and so reduces 
the g-g lensing signal. Despite this effect, we are able to re- 
produce the fraction of satellites in the overall sample rather 
well: as shown in tabled] the best-fit satellite fractions are 
in most cases within 10% of the correct value. This quantity 
therefore appears to be relatively robustly extracted from 
the data. This works best if we simply ignore the h-h term 
and only fit out to 2h~ Mpc. In a more general fit with the 
h-h term and out to larger radii, we find there is some degen- 
eracy between the components, and the extracted satellite 
fraction can be too low in some cases (table |3J. 

For the halo model, we have to choose the radial distri- 
bution of galaxies, which we model as a NFW profile with 
its own concentration parameter c g . The concentration of 
subhalo radial distrib ution depends sensi t ively on how the 
subhalos are selected. iNagai fe Kravtsovl <l2005l) show that 
selection of subhalos by their maximum circular velocity 
Vmax results in concentrations of subhalo number density 
profiles in clusters in the range ~ 1.5 — 4, while selection us- 
ing luminosity or stellar mass results in considerably larger 



concentrations (c ~ 3 — 10). Although luminosities have 
been assigned to the subhalos in our catalogs, the assign- 
ment did not take into account the position of the subhalo 
within its host and depended only on the subhalo's Vmax. 
We can therefore expect that the radial distribution of the 
cluster subhalos in simulations will be characterized by con- 
centrations similar to those measured for subhalos selected 
using Vmax (i.e., c g ~ 1.5 — 4). Results of our analysis over 
a broad range of host halo masses, between 10 12 /i _1 Mq to 
10 14 /i _1 Mq, indeed suggest values of c g significantly lower 
than Cd m - Table 13 gives the true c g value from simulations 
weighted over the halo mass distribution for a given lumi- 
nosity bin. 

It should be noted that shallow profiles for the satel- 
lite distribution are not a unique prediction of the simu- 
lations. Semi-analytic simulations find c g ~ Cd m by identi- 
fying galaxies before merging into larger halos and follow- 
ing their identity even after their surro unding dark matter 
has been stripped off llGao et alj l2004l. Observations sug- 
gest that the galaxy distribution can be well described by 
the NFW profi l e in rich clusters (e. g . , ICarlberg et alj Il997l : 
iLin et aljl20ol lHansen et al.ll2004l) These works also find 
that galaxies have concentrations significantly lower than 
expected for their parent halos. Since the actual value of 
c g depends on the physical processes that affect the specific 
type of galaxy under consideration, we treat c g as a free 
parameter to be determined by fitting the simulations (or 
data) . 

We find that in order to fit the simulations, we need 
the satellite distribution to be shallow (table 0, consistent 
with the value found by directly fitting to the satellite dis- 
tribution in the same simulations, except for the two high- 
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Table 3. Satellite fraction and concentration parameter of satellite radial distribution from simulations and from the halo model fits for 
fits out to 2h~ 1 Mpc without h-h term and out to 13/i -1 Mpc including the h-h term. 



Luminosity a tI ue a flti2h -i Mpc a flt ,i3h-iM P c c 9,tru C c g flt 2h -i Mpc c Sifitjl3h _i Mpc Scatter in V max - L 



L3 


[-19,-20] 


0.23 


0.26 ± 0.03 


0.22 ±0.02 


1.5 


2.2 ±1.1 


3.2 ±0.4 


no 


LI 


[-20,-21] 


0.20 


0.21 ±0.01 


0.18 ±0.02 


3.0 


4.7 ± 1.8 


5.9 ±0.7 


no 


L5 


[-21,-22] 


0.15 


0.15 ±0.02 


0.10 ±0.02 


3.5 


1.5 ±1.3 


2.1 ±0.4 


no 


L3 


[-19,-20] 


0.23 


0.26 ±0.03 


0.22 ±0.02 


1.5 


2.8 ±1.4 


3.9 ±0.5 


yes 


LI 


[-20,-21] 


0.21 


0.20 ±0.02 


0.17 ±0.02 


2.5 


10 ±4 


12 ±1 


yes 


L5 


[-21,-22] 


0.19 


0.16 ±0.05 


0.16 ±0.03 


2.3 


20 


20 


yes 



est luminosity bins. As discussed above, the low value of c g 
found in simulations does not imply that the same value has 
to be appropriate to fit the data: the radial distribution of 
galaxies of a given luminosity inside larger halos is affected 
by star formation and dynamical processes both before and 
after falling into the larger halos, which is likely to be quite 
complicated. What is important for the current discussion 
is that the halo model should reliably determine the radial 
distribution of galaxies from the g-g signal, allowing one to 
determine it from observations. 

Table El also shows the best-fit value of c g from the full 
signal. As shown, even with the ten times reduced SDSS 
errorbars, we are not highly sensitive to c g . The result for 
L5 with scatter is quoted without errorbars because the fit 
program chose the maximum value allowed in the fits, with 
very small errorbars. This result (and the high value for L4 
with scatter) may be an indication that AE averaged over 
such a large range of halo masses does not lend itself well 
to this type of fitting for a single value of Mo and c g . Once 
again, this result emphasizes that the fitting is most mean- 
ingful for groups of lens galaxies chosen to have a narrow 
distribution in central halo mass, which may not necessarily 
be the case for r-band luminosities even with narrow (1- 
magnitude wide) bins. The results of the fits suggest that 
the radial distribution can in principle be extracted from 
the observations, although we note that when applying this 
to the actual data, the observational errors were found to be 
large and unable to provide a strong co nstraint from the cur- 
rent data samples fS eliak et al.ll2005h . Note that when the 
fit was performed to 13/i _1 Mpc using the combined signal, 
the bias values b(Mo) were found to be lower than the values 
in table |21 and an examination of figure Q suggests that on 
intermediate scales (2-4 /i _1 Mpc), there is some interplay 
between c g and fe(Mo) that is causing the high best-fit val- 
ues of c g and the low values of b(Mo). A better approach 
may be to determine the bias b(Mo) using the central signal 
alone as in table|5| then fix it to that value in the fits for the 
full signal rather than allowing it to vary; this will encourage 
lower, more reasonable values of c g . 

We note here that the satellite radial distribution is 
more important for the overall shape of the satellite contri- 
bution to g-g lensing signal than the halo occupation num- 
ber as a function of mass (N(M)). Changes in the slope e in 
(N(M)) oc Af £ mostly change the amplitude of the satellite 
signal without changing its shape (see figure 7 in GS02). This 
result is expected, since the satellite signal is dominated by 
halos with masses in the range (10 13 — 10 14 )/i _1 Mq . Change 
in the slope changes their relative abundance, but more or 
less preserves the radial shape of their signal. Consequently, 



when we attempted fits for both e and a, the slope e was 
almost completely degenerate with the fraction of satellites 
a, but the actual halo mass probability distribution in this 
mass range was less affected. Hence, as discussed already, we 
did not fit for e, but instead assumed fixed e = 1 throughout. 

Another aspect of the halo model that can be tested 
with simulations is the amount of dark matter attached to 
subhalos within the larger halo, though as discussed, we only 
approach this problem in an average sense. For g-g lensing, 
the tidal stripping of the outer layers of dark matter attached 
to the satellite is of little importance, since the lensing sig- 
nal of that component is swamped by the host halo dark 
matter signal. In GS02 it was simply assumed that all of 
the dark matter mass is attached to the subhalo, so that in 
the inner region, the lensing signal of the satellites is the 
same as of the central galaxies of the same luminosity. Fig- 
ure which shows the noncentral signal for one luminosity 
bin, with the best-fit signal assuming unstripped, moder- 
ately stripped, and completely stripped subhalos, shows that 
this is a good approximation in the inner parts. One should 
perhaps not be too surprised by this agreement, since by 
construction in these dissipationless simulations, the galax- 
ies sit on top of subhalos that have survived all the merging 
and tidal stripping inside the halos. Reality could be more 
complicated: one could have dark matter in subhalos en- 
tirely stripped while the more compact stellar component is 
preserved. 

While the agreement between the simulations and halo 
model is already good, one can improve the agreement fur- 
ther by modeling the transition between the satellite lensing 
signal and the host lensing signal, which occurs around 100- 
200 /i _1 kpc. We assume that only the dark matter within 
0.4 of original virial radius is attached to the satellite (fig- 
ure|nj. This requires that on average 50% of the dark matter 
is stripped from the satellites. We find the signal to be rela- 
tively insensitive to the exact value of truncation radius and 
consequently g-g lensing cannot be used as a strong probe 
of tidal stripping process; however, as shown in figure |S| 
we can exclude the extreme possibility that all halos are 
fully stripped. To study this issue observationally, it is best 
to select a galaxy sample with a particularly high satellite 
fraction a, such as the galaxies in dense enviroments. 



5 CONCLUSIONS 

The goal of this paper is to present a detailed comparison of 
halo model and simulation predictions of galaxy-mass corre- 
lations as observed in weak lensing. In its simplest form, and 
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L4, without scatter, satellite contribution 
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Figure 6. AS(r) for satellite contribution to AS from L4 with- 
out scatter. The points are from simulations, and the lines are 
(as labeled) the best-fit signal with no mass stripped, partially 
stripped (truncation at 0.4rj,i r ), fully stripped, and a mixture of 
half partially stripped and half fully stripped. 

for a given cosmological model, the halo model has 3 free pa- 
rameters, which depend on the physics of galaxy formation. 
These are the typical halo mass of central galaxies (where 
"typical" is between the median and mean halo mass), the 
fraction of galaxies that are satellites, and the radial distri- 
bution of satellites inside the larger halo. 

The halo model provides a fairly accurate description 
of the g-g lensing signal from simulations, in the sense that 
it can determine the virial mass distribution of central and 
non-central galaxies and their relative ratios, as well as the 
radial distribution of both galaxies and dark matter. One 
must be careful when interpreting the halo masses from 
NFW fits if the halo mass distribution is broad as is ex- 
pected for luminous galaxies (L > few L*). In this case, g-g 
lensing is determining something between the mean and the 
median mass and may not recover the correct radial distri- 
bution of the satellite galaxies in their host halos. If one is 
interested in the mean halo mass, as in the case of appli- 
cation to bias studies JSeliak et al.ll200fjl . then an upward 
correction needs to be applied; this correction is small at the 
faint end but can be quite significant at the bright end, es- 
pecially if there is a significant scatter in the luminosity-halo 
mass relation. 

There are additional parameters that can be added to 
the description, such as the generalization of the assumed 
halo mass probability distribution, but they do not improve 
the fit and are strongly correlated with these 3 parame- 
ters. Galaxy-galaxy lensing also allows a determination of 
the halo mass profile of both central galaxies and satellites. 
In the latter case, one can constrain the amount of satellite 
mass stripping inside larger halos, which determined how 
much dark matter remains attached to the satellites. To in- 
vestigate the radial profiles of dark matter it is best if one 



identifies samples that consist predominantly of central or 
satellite galaxies, respectively. 

Galaxy-galaxy lensing is also quite sensitive to the dark 
matter concentration parameter Cd m , which (in standard 
cold dark matter models) is not a free parameter but is fixed 
for any given cosmological model. It varies modestly as a 
function of cosmology, so in principle one could use this fea- 
ture to determine the cosmological model. For example, for 
a typical L* galaxy with mass around (1-2) x 10 12 /i _1 Mq, 
we have Cd m = 12 in a ag = 0.9, fl m = 0.3 cosmology, which 
drops to Cdm = 10 if as — 0.7 is used instead. This change is 
small and we find that it is degenerate with the amount of 
scatter, which also leads to a lower value of concentration pa- 
rameter. With current samples, one cannot yet reliably dis- 
tinguish between these models, but one can test the overall 
consistency. The accuracy and reliability with which the con- 
centration can be extracted from the data may be improved 
if galaxies in underdense regions, for which the satellite frac- 
tion is lower, are selected. In this case, one can use the signal 
to a larger distance since it is not swamped by the noncentral 
component. Concentration can also be determined more ac- 
curately with deeper survey s, which can probe d ark matter 
halo at smaller separations IIHoekstra et al]|2004T) . although 
modeling angular projection effects in the absence of reliable 
redshift information for lens galaxies may be difficult. 

We find that the data are rather insensitive to the trun- 
cation radius of the satellites inside the larger halo as long 
as it is not very small. We can, however, test the extreme 
possibility that satellites have no dark matter attached to 
them. This is unlikely to be the case for all satellites even 
though it may be v alid for some frac tion of those near the 
center. For example. iGao et al] i2004) argue that up to 40% 
of galaxies in clusters may not have associated DM halos, 
a conjecture which may be testable with the future weak 
lensing measurements (see Fig. |SJ . 

Additional information about the galaxy-dark matter 
halo connection may be more challenging to extract from 
the data. For example, we have assumed that the number of 
satellites inside the halo scales linearly with halo mass. If this 
assumption is dropped and a more general power law rela- 
tion is allowed, a degeneracy is developed between the satel- 
lite fraction and the power law slope, but the overall fraction 
of galaxies in (10 13 — 10 14 )/i _1 Mq halos is preserved. 

We note here that the power law slope and other pa- 
rameters may be constrained by other observations such as 
the galaxy auto-correlation function analysis, velocity infor- 
mation in redshift surveys and direct counting of galaxies in 
groups and clusters. In fact, a halo model similar to the one 
use d here was applied t o the galaxy auto-correlation analy- 
sis iZehavi. et al.ll2004Tl . Combining and comparing the two 
data sets should provide important consistency checks and 
will further constrain the galaxy-dark matter connection. 
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